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Abstract 

We consider the raise and peel model of a one-dimensional fluctuating interface in 
the presence of an attractive wall. The model can also describe a pair annihilation 
process in a disordered unquenched media with a source at one end of the system. For 
the stationary states, several density profiles are studied using Monte Carlo simulations. 
We point out a deep connection between some profiles seen in the presence of the wall 
and in its absence. Our results are discussed in the context of conformal invariance 
(c = theory). We discover some unexpected values for the critical exponents, which 
are obtained using combinatorial methods. We have solved known (Pascal's hexagon) 
and new (split-hexagon) bilinear recurrence relations. The solutions of these equations 
are interesting on their own since they give information on certain classes of alternating 
sign matrices. 
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1 Introduction 



The raise and peel model (RPM) is discussed in several papers [U [2] (for a review see 
[3]). This is a one-parameter dependent stochastic model of a one- dimensional fluctuating 
interface. The RPM can be seen as a model for wetting in which the adsorption is local (the 
interface raises) and the desorption is nonlocal (the interface is peeled). The parameter is 
given by the ratio of the adsorption and desorption rates. When the two rates are equal (the 
Razumov-Stroganov point) the Hamiltonian which describes the time evolution of the system 
is conformal invariant and the stationary states have remarkable combinatorial properties. 
The origin of these properties can be traced back to the definition of the model in terms 
of generators of the Temperley-Lieb (TL) algebra at the semigroup point. This algebra has 
two representations which are relevant for us: one in terms of link patterns (which makes 
the connection with the interface model) and another one in terms of an integrable XXZ 
quantum chain with L sites (which allows one to obtain the spectrum of the Hamiltonian 
and its conformal properties). 

The fluctuating interface is described by RSOS (Dyck) paths which can be interpreted 
as clusters of tiles deposited on a substrate. The heights at both ends of the system vanish. 
The clusters touch each other at contact points. In [1], a conjecture for the average number 
of clusters at the Razumov-Stroganov point in the stationary states was given. According to 
this conjecture, the average size of a cluster increases like L^/^ (L is the size of the system, L 
even). In the present paper we give combinatorial arguments in favour of this conjecture. We 
also show that the average size of the cluster at one of the ends of the system also increases 
like L^/3_ r^Yiis is in contrast to a model in which the RSOS paths are taken with equal 
probabilities when one obtains L^/^. 

An interesting result obtained in [2] using Monte Carlo simulations is the behavior of the 
density of contact points in the finite-size scaling limit. It has the functional form expected 
from conformal invariance for a one-point function with an unexpected exponent. Further 
details on this observation will be given in the next sections. 

In the present paper we extend the RPM at the Razumov-Stroganov point in two ways. 
Firstly we consider the case when the size of the system L is an odd number. This implies 
a change of the configuration space of the model, one has ballot paths with fixed ends: one 
at height 1 the other one at height zero. We call them one-step Dyck paths. The model 
is still based on the TL algebra at the semigroup point. We also consider the RPM in the 
presence of a wall (RPMW). The model is based on the one-boundary TL algebra |5l [6] at 
the semigroup point. The configuration space is now given by ballot paths [7] with one end 
fixed at zero and the other one free, corresponding to the wall. In the bulk the rules for 
adsorption and desorption as in the original RPM but at the wall one has a new adsorption 
process with an arbitrary fixed rate a (for convenience the bulk rates are taken equal to one). 

For both extensions conformal invariance is obeyed and the stationary states have new 
combinatorial properties. 

There is an alternative interpretation [8] of the RPMW. One can see the stochastic 
process as a pair annihilation processes in an unquenched disordered media (see Appendix 
A for detailed explanations). In this description, one has clusters of tiles which either touch 
each other or are separated by impurities (defects). There are no empty sites. The defects 
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can either hop over a cluster peehng its surface or they can annihilate in pairs if the defects 
are on neighboring sites. The clusters change in shape and number like in the usual RPM. 
The rates of all these processes are fixed at the value one. At one end of the system one has 
a source of defects which acts with a rate a. The source acts as follows: if a cluster touches 
the end of the system, a layer of the cluster is desorbed and a defect is added at each of 
the two ends of the peeled cluster. In the stationary states, the average number of defects 
corresponds to the average height of the ballot paths at the wall. 

In the " defects" interpretation of the model, the RPM for L odd can be seen as a system 
in which one impurity (defect) hops at large distances in an unquenched disordered media. 
During the hop, the media is changed (the cluster over which the defects hops, is peeled). 

This paper has two distinctive parts: a physical one and a mathematical one. Physics is 
in the first part of the paper (Sections 2-4), mathematics is concentrated in the second part. 
The two parts can be read independently. 

In Section 2 we present the stochastic models describing the RPM models with or without 
a wall. This presentation is suitable for Monte Carlo simulations. We also define some 
relevant observables. The alternative description of the model in terms of defects is given in 
Appendix A. 

In Section 3 we present without proof some exact results for the stationary states obtained 
in the last part of the paper. We give the expressions for the average number of clusters 
for any size of the system. In particular we give the large L behaviour. These expressions 
are relevant because they are the integrated quantities of the density of contact points, the 
latter being obtained from Monte Carlo simulations. In particular we derive the values of 
the critical exponents. 

We show a surprising identity between the probability density function (PDF) to have 
the first cluster at distance x from the wall for the RPMW and the PDF to have defect at 
distance x from the boundary for the RPM for an odd number of sites. As shown in Section 
4 this PDF has remarkable properties in the finite-size scaling limit. 

We also give some conjectures for the average values of the number of sites where adsorp- 
tion can take place. From their values one can derive the average number of tiles desorbed 
(desorption takes place through avalanches [2]). 

In Section 4 we present the results obtained from Monte Carlo simulations. Firstly we 
remind the reader about the observation [2] that in the stationary state of the RPM (L even) 
the density of contact points, in the finite-size scaling limit is given by the expression expected 
from conformal invariance with an exponent 1/3 while we should expect an exponent 2/3 
(the one-point function of an operator with conformal spin vanishes). 

We next consider the presence of the wall in the model with a boundary rate a and study 
the density of contact points profile in the finite-size scaling limit. It is independent of the 
boundary rate a and has a functional form suggested by conformal invariance with the same 
exponent 1/3 as in the absence of the wall. 

We also study the density of defects profile and find that in the finite-size scaling limit 
it has the functional dependence expected from conformal invariance with (again!) an un- 
expected exponent. 

In the finite-size scaling limit, the probability distribution function to have a cluster at 
distance x/L from the wall stays, surprisingly, a probability density function (its integrated 
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value is equal to one) and is independent of a i.e. is universal. The functional dependence of 
this PDF is similar to what is expected from conformal invariance for a one-point function. 
No explanation for this observation was found. 

In Section 5 we present the derivation of the RPM and RPMW in continuum time 
using the Temperley-Lieb algebra [9] and its one-boundary extension [5] at the semigroup 
points and make the connection with the spin one-half XXZ quantum chain. Based on this 
connection and the integrability of the chain, one can show that our systems are conformal 
invariant. This section is a review of known results. 

In Section 6 we present firstly the bilinear recurrence relations called the Pascal's hexagon. 
The relations are not new but the solutions, specified by the boundary conditions are. We 
next propose some new bilinear relations that we called the split-hexagon relations and give 
their solutions for certain boundary conditions. This section is pure mathematics. 

In Sections 7 and 8 we make the contact with physics and derive the results enumerated 
in Section 3. Firstly in Section 7, for the RPM we derive the probabilities to have k clusters 
for a system of size L {L even and odd). 

In Section 8 we consider the RPMW. For L odd we show a remarkable connection be- 
tween the RPMW and the RPM for the same value of L. In particular one can show that 
independently of a value of the boundary rate a all properties of clusters, except the first, 
are identical. The case of the RPMW with L even and arbitrary values of a is more subtle. 
Nevertheless for the case a = 1 we were able to derive the probability of having k clusters 
in a system of size L. 

Our conclusions are presented in Section 9. 



2 Raise and peel models with and without a wall 

We consider a one-dimensional lattice with (L + l) sites. An interface is formed by attaching 
at each site non-negative integer heights hi which obey the restricted solid-on-solid (RSOS) 
rules: 

h,+i-hi = ±l, hi>0, 2 = 0,1,...,L. (2.1) 

We will consider three kind of interfaces, each with different configuration spaces depending 
on the conditions at the boundaries (the values of ho and h^): 

a) Dyck path configurations. 

One takes L even and ho = hj^ = 0. There are 

L\ 

^ (L/2)!(L/2 + 1)! ^ ■ ' 

configurations of this kind. An example of such a configuration for L = 16 is shown in 
Figure 1 . 

To characterize the interface, it is useful to define several quantities. A contact point is 
a site j where hj = 0. A cluster is the domain between two consecutive contact points. 
The size of the leftmost cluster is j, if j 7^ is the smallest number for which hj = 0. 
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Figure 1: (Color online) A configuration of the interface given by a Dyck path with three 
contact points and two clusters for a lattice size L = 16. The size of the leftmost cluster 
is equal to 8. The substrate is also shown (the dashed region). 



1 7 L=15 

Figure 2: (Color online) A one-step Dyck path configuration with four contact points 
and four clusters for L = 15. The leftmost cluster has size 7. 




It is useful to visualize the interface as a film of tiles (tilted squares) deposited on a 
substrate defined by the Dyck path hi = ^~'^2^^ (z = 0, 1, . . . , L). In this picture the 
clusters can be seen as droplets of a fluid deposited on the substrate. 

b) One-step Dyck paths configurations. 

One takes L odd, = 1 and hi = 0. There are Cl+i configurations of this kind. An 
example of such a configuration for L = 15 is shown in Figure 2. 

The one-step Dyck paths can be mapped onto configurations with clusters and one 
defect (impurity). In order to obtain these configurations one draws a horizontal line 
at height one through the first (leftmost) cluster starting at ho = 1. If the first cluster 
has size x {h^ = 0), the horizontal line intersects the cluster at (x — 1). One puts a 
defect at the point (x — 1/2) and one lowers the cluster by one unit: hj {hj — 1) 
Vj = 0, 1, ... ,x — 1. This procedure is illustrated in Figure 3 for the configuration 
shown in Figure 2. 

c) Ballot paths. 

One takes (for both L even and odd) h^ = and /iq free {ho = L, L — 2, . . . , or 1) 
There are 

L \ 



(2.3) 



configurations of this kind (here \n\ is the integer part of n). An example of such a 
configuration for L = 14 is shown in Figure 4. 

The ballot paths configurations can be seen as droplets of a fluid in the presence of a 
wet wall. 



^ In the usual definition , a ballot path is defined by taking Hq fixed. 



5 




1 7 L=15 ^ 1 6 7 L=15 



Figure 3: (Color online) Mapping the one-step Dyck path shown in Figur^S] onto a 
configuration with one defect and clusters. 




1 6 11 L=14 



Figure 4: (Color online) A ballot path for L = 14. One has one contact point and one 
cluster of size 14. 



Similar to the case of one-step Dyck paths, one can map the ballot paths onto config- 
urations with clusters and defects (impurities). We will call them configurations with 
defects. The mapping follows the same procedure as the one used for one-step Dyck 
paths and is illustrated in Figure 5. Notice that the height at the origin ho is equal to 
the number of defects. 

The configurations with defects can also be seen as an interface in which the droplets 
have empty spaces between them. 

In the raise and peel model (RPM) the dynamics of the interface is described in a trans- 
parent way in the language of tiles (tilted squares) which cover the area between the interface 
and the substrate. We consider the interface separating a film of tiles deposited on the sub- 
strate from a rarefied gas of tiles. The interface can be Dyck path, one-step Dyck path or 
ballot path configurations (the dynamics for defect configurations is discussed in Appendix 
A). 

The evolution of the system in discrete time (Monte Carlo steps) is given by the following 
rules. With a probability Pj = 1/(L — 1) a tile from the gas hits the site i, i = 1, . . . , L — 1. 
Depending on the value of the slope Si = (/ij+i — /ij_i)/2 at the site i, the following processes 
can occur: 

i) Si = and hi > hi-i. 

The tile hits a local peak and is reflected. 

ii) Si = and hi < hi^i. 

The tile hits a local minimum. With a probability 1 the tile is adsorbed [hi hi + 2). 

iii) Si = l. 
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Figure 5: (Color online) The mapping of the ballot path shown in Figure 4 onto a config- 
uration with four defects. 



o 




12 8 i=16 



Figure 6: (Color online) A desorption event. The incoming tile at site 1 triggers 
an avalanche of 5 tiles, which are shaded. All the shaded tiles are removed in the 
desorption event. 



With probability 1 the tile is reflected after triggering the desorption of a layer of tiles 
from the segment {hj > hi = hi^j,, j = i + 1, . . . ,i + b — 1), i.e. hj hj — 2 for 
j = i + 1, . . . ,i + b — 1. This layer contains b — 1 tiles (this is always an odd number). 
For an example see Figure 6. 

iv) Si = -1. 

With probability 1 the tile is reflected after triggering the desorption of a layer of tiles 
belonging to the segment {hj > hi = hi^b, j = i — b + 1, . . . ,i — 1), i.e. hj hj — 2 
for J = i — 6 + 1, i — 1. 

Notice that the adsorption and desorption rates were taken equal, this is the RPM at 
the Razumov-Stroganov point. The model was studied in detail in the case of Dyck path 
configurations (L even) also when the adsorption and desorption rates are different and the 
phase diagram of the model was obtained (see Refs.[ll[2], for a review see |3]). The Razumov- 
Stroganov point is special in two ways. Firstly, the Hamiltonian which gives the continuous 
time evolution of the system can be mapped into an XXZ spin one-half quantum chain which 
is integrable [10] . The finite-size scaling limit of the Hamiltonian spectrum is given by c = 
Virasoro characters (c is the central charge of the Virasoro algebra) and hence the system 
is conformal invariant. Secondly, the probability distribution function (PDF) describing the 
stationary state of the system for finite systems has remarkable combinatorial properties 
[m [12] which allows one to obtain exact results for physical observables. 

The RPM was also considered in the L odd case and it was shown that the defect (see 
Figures) makes Levy fiights and behaves like a "relativistic" random walker ( dispersion 
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Figure 7: (Color online) The adsorption of a half-tile at the first site and desorption of a 
layer touching the boundary for the ballot path shown in Figure 4. 

relation (x^) ~ t^) [8]. In the present paper we are going to present more results which will 
show some surprising properties of this model. 

The PDF describing the stationary states is expressed in terms of Dyck paths only. 
The dynamics seen in the space of defect configurations is interesting: it describes the pair 
annihilation of two defects in an unquenched random media. The average density of defects 
decreases like 1/t |8]. 

We are going to extend the RPM model acting on ballot paths by making the first site 
active (site in Figures 4 and 5). We will define in this way the RPM in the presence of a 
wall (RPMW). 

The evolution of the system (Monte Carlo steps) is given by the following rules. With 
a probability Pi = 1/(L + a — 1) a tile from the gas hits the site i, i = l,...,L— 1. The 
changes of the interface produced by the hits are the same as in the RPM model. With a 
probability Pq = a/{L + a — 1), a half-tile hits the site 0. The boundary rate is equal to a, 
as opposed to the bulk rates which are equal to 1. If the slope Sq = hi — is equal to 1, the 
half-tile gets adsorbed (see Figure?). If Sq = —1, the half-tile is reflected. Notice that the 
adsorption process is local. This is not the case if one uses the defect configurations picture 
(see Appendix A). 

The RPMW in its time continuous version, can be obtained in a simple way using a 
Hamiltonian expressed in terms of generators of the one-boundary Temperley-Lieb algebra. 
This is explained in Section 5. Similar to the RPM model, the stationary states of the RPMW 
have magical combinatorial properties for finite values of L [T3l [TU [15] . New expressions 
and relations are summarized in Section 3 and derived in Sections 7 and 8. Moreover, the 
Hamiltonian can be mapped onto an integrable XXZ spin one-half quantum chain and the 
finite-size scaling limit of the Hamiltonian spectrum is known [TB]. It can be expressed 
in terms of Virasoro characters therefore, similar to the RPM, the RPMW is conformal 
invariant. The consequences of this observation will be presented in Section 4. 

3 Some exact results and some conjectures 

In this section we present some exact results which will be discussed in Sections 7 and 8. 
These results besides being interesting on their own, will be used in Section 4 to derive 



8 



critical exponents and to check the data obtained in Monte Carlo simulations. We consider 
the RPM with and without a wall at the Razumov-Stroganov point and the boundary rate 
a. 

We first introduce two useful functions S{L,n) and P{L,n,m). They are defined as 
follows: 



S(L,„) = 2-'""Yl7:rr^ U 



L-3-J f T _ \ n+p 



(L- m -p-l)\ 



\{2p-l)\\ ^J-^ {L-2n + 3p)\ 



X 



L 3 ^ (2L + 2n-6p-3)!! 



n (2z,_2[!^J+4p+l)!!' ^^-^^ 



\ ofT i\ (2n - m)!(n - 1)! 
P(L,n,m) = 21- 2 J 5'(/, 77, _ 1) LA 

^ ' ' ^ ^ ' ^ m\{n - m)\{2n - 1)\ 

{L-2n + m- 1)!(L - n)\{2L - 2n + 2m - 1)!! 



(3.2) 



(L-2n-l)!(L-n + [f J)!(2L - 2n + 2[^J - 1)!! ' 

The expressions (13. ip have first appeared in [131 E] • The expressions for 

S{L) = S{L, [{L-l)/2\) 

appear also in the problem of enumeration of vertically symmetric and vertically and 
horizontally symmetric A^^ pxp alternating sign matrices (a known topic in combinatorics 
[m [TBI [19]): 

^ ' ^""^^ o<iLi (4^ + 2)!(4z + 3)! ^ ' 

S{2N — 1) = A^^_^/A2N-i = ^AN+i/^2N+i = II + 1)! ' ^'^'^^ 

We next fix the notations to be used below: for the RPMW of a size L and boundary rate 
a an observable $(. . .) is denoted as $i^''(. . .). For the RPM of a size L the same function is 
denoted as $l(. . .). Notice the obvious relation $^°''(. ..) = $/,(...). 

We now state the following results for the stationary states of the RPM and RPMW. 

In the RPM, the probability of having k clusters for a system of size L is: 

PAh) = -^P(L,L^J,A;-l). (3.5) 

For L even this result coincides with a conjecture made by de Gier [1]. 

In the presence of a wall (RPMW), the probabilities to have k clusters for a system of 
size L are: 

atutty: ^i'^^) ^ ^-^^) ^ ^P(L,^,.-1), (3.6) 
^even, ^(i)^^^ ^ ^{p(L + 1, f , A.) + P(L + 1, f - 1, A: - 1)}. (3.7) 
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Note the remarkable fact that for L odd, the probabihties to have k clusters for the RPM 
with or without a wall are the same. This statement is valid for any value of a. Another 
connection between these two systems (with or without wall for L odd) will be presented 
below. The relation (13. 7p for L even is valid only for a = 1. 

Knowing the probabilities to find k clusters, one can compute the average number of 
clusters (/c)^ and {k)j^\ for the RPM and RPMW, respectively. We give here only the large 
L limit behavior of these numbers: 

r r(iZ3)V3L2/3^ 0.738488 L2/3, (L even) 
{k)L ^-^ - (3.8) 

[ aL2/3, (L odd) 

(^)l'' — ^ aL'^^^, (L even or odd, a arbitrary) (3.9) 

where approximants give for a the value a ~ 0.5056(3). Observe that in all cases the average 
size of a cluster increases like L^^'^. Notice that for the RPM for L even one has the largest 
number of clusters. The explanation is very simple. In the case L odd or in the presence 
of the wall, the first cluster near the active boundary has a size of the order of L (see the 
comment after eg. (13. 151) ) which leaves less available space for the other clusters. 

We now consider the probability densities Fj^^\x) to have the leftmost (first) cluster end 
at a distance x from the wall for both L even and odd. These probability densities are 
related to the probability density Dl{x) to have the leftmost cluster end (which coincides 
with the position of the defect) at the point x for the RPM for L odd: 



L odd, 
a arbitrary: 



Fi''\x) = Dl{x), x = 1,3,...,L, (3.10) 



F«(x) = Dz.+i(x + l), x = 0,2,...,L. (3.11) 

Moreover, the conditional probabilities to have the first cluster at distance x and any given 
configuration to the right of it are the same in the two models for any value of a if L is odd, 
and for a = 1 if L is even. In particular, the densities of the contact points (definded on 
page m and illustrated in Figures 1, 2 and 4) at a distance x in the RPMW and in the RPM, 
Nj^\x) and Nl{x), respectively, satisfy the relations 

aa*itrty: ^^i^'W = A'.W. .= 1.3....,L. (3.12) 

Leven, Ni'\x) ^ N:,^,(x + 1), a: = 0. 2, L. (3.13) 

(Jj — -L . 

As discussed in Appendix A, for L odd, in the RPM to a one-step Dyck path corresponds 
the one defect picture in the configuration space. The probability density to have the defect 
at a distance x is obviously the same as at (L — x), a; = |, |, . . . , (L — i). This observation 
implies the following symmetry relations: 

o atbtty: ^ fI'\l ^ 1 - .) , .= 1,3,....L. (3.14) 
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Fi'\x) = Fi'\L-x), x = 0,2,...,L. (3.15) 

The symmetry relations fl3.14p and flS.lSp are surprising since one could expect that the 
probability density function should be biased towards the wall. Another important conse- 
quence is that the average size of the first cluster is L/2. This is in contrast with the RPM 
for L even when it is of order L^/^ (see Section 4). 

Before closing this section we give a conjecture for the RPMW with the boundary rate 
a = 1. The fraction of the interface where adsorption can take place (this is the average 
number of local minima of the interface divided by L) is 



4" 



6L2 + 8L - 5 
4(2L + l)(2L + 3)' 

(6L2 + 14L + 9)(L-1) 
4L(2L + l)(2L + 3) 



(L even) 

(3.16) 

(L odd) 



These expressions were checked up to L = 9. It is amusing to note that simple guesses are 
possible for A similar situation occurred also for the RPM [1]. The expressions (13.161) 
will be used in Section 4 to characterize the avalanches occurring in the model. 



4 Density profiles in the raise and peel models with 
and without a wall in the stationary states 

As mentioned in the previous sections, in the continuum limit, the RPM and the RPMW 
are conformal invariant and are described by a c = conformal field theory [lOl |16] . If the 
central charge of the Virasoro algebra vanishes, the scaling indices (highest weights of the 
irreducible representations of the Virasoro algebra) are: 

A., ^ <^^^^|f^. (4.1) 

where p and q are nonnegative integers. This implies that Ap^q is either an integer or it is 
equal to 1/3 plus an integer. 

A local operator </> is characterized by its scaling dimension X and conformal spin s: 

X = Ap^q + Api^qi, S = Ap^q — Api ^qi . (4.2) 

If X is the distance from the origin (which coincides with i in the discrete versions de- 
scribed in Section 2, i = corresponds to x = 0) and L is the size of the system, in the scaling 
limit, one expects [20] the following density profile in the case of symmetric boundaries "r": 



where C is a constant. 



(0(x/L)),,, = 7^-^^^, (4.3) 
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Figure 8: Typical configuration in the stationary state for the RPMW with a boundary rate a = 1. 
The system has a size L = 128. hi is the height at site i. 



The profile of a primary operator vanishes if its conformal spin is not zero [20] therefore 
one should have 

X = 2Ap,,. (4.4) 

The fact that one can make predictions for the functional dependence of the density 
profiles in stationary states is a consequence of the conformal invariance of the Hamiltonian 
which gives the time evolution of the stochastic process. 

If the two boundaries are different (r 7^ s), one expects |20j : 



i<P(-/L))r. = ^;"^7i4l\ (4.5) 
[Lsin(-)] 

where the function $r,s is calculable. 

Before discussing density profiles, it is instructive to have under our eyes typical profiles 
for finite lattices. For the RPM, L even, typical profiles are shown in Ref. [2]. In the case of 
the RPMW and finite L, the profiles depend on the boundary rate a. In Figure 8 we show 
for a = 1, a typical configuration for L = 128. One notices that ho = 4: which implies that 
in the corresponding defect configuration one has 4 defects (see Appendix A). One has a 
large leftmost cluster of size 124. In this paper we will not discuss the profile of the heights 
(for the RPM (L even) this was done in full detail in [2]) but we will show that the average 
height at the origin behaves like (ho) ~ ln(L). 

We discuss several density profiles. 

1) Density of contact points. 

a) The RPM model (L even). Density of contact points in Dyck path configurations. 

Since the average number of contact points is known (see (13. 8p ). if (14. 3 p is valid, we have 
to choose X = 1/3. The functional dependence of the density is therefore also fixed as well 
as the constant C in (14. 3p : 

C= - ^r(4) ^ 0.753149. (4.6) 

In order to check the prediction of conformal invariance, using Monte Carlo simulations 
on large lattices, we have measured the density of contact points Nl{x). In Figure 9 we 
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Figure 9: (Color online) The RPM for L even. The density of contact points Nl{x) 
multiplied by [Lsin(7rx/L)/7r] 1/^0-753149 as a function of x/L for L = 512, 1024, 2048 
and 4096. 



show A'"i(x)[Lsin(7rx/L)]i/'^/C for several lattice sizes. If the prediction (14.31) is correct, one 
should obtain the value 1 for this quantity. This is indeed the case. (Data with poorer 
statistics but indicating the same result were shown already in [2]). 

One does not need to find the density of contact points for the RPM in the case L odd 
since it coincides (see Section 3) with the one observed in the RPMW for a boundary rate 
a = 1 which is going to be discussed next. 

b) The RPMW model. Density of contact points in ballot path configurations. 

In this case the boundary conditions at x = and x = L do not coincide since for 
X = one has the wall and for x = L, one has hi = 0, fixed. For this reason one expects an 
expression like the one given by (14. 5p [20] with the function ^w,o{cos{nx / L)) to be determined 
(here w corresponds to the wall, corresponds to Hl = 0). 

We start with the case in which the boundary rate is a = 1. Since the average number 
of clusters is known exactly (see (13.91) ). the exponent X in (14.31) is X = 1/3. Monte Carlo 
simulations on large lattices suggest the following simple expressions for $^ o(cos(7rx/L)): 

$^,0 = /?[l-cos(f)]^/^ (4.7) 

where 

(3 = ^/3a ~ 0.8757(8). (4.8) 

This number is compatible with the data shown in Figure 10. 

In order to make sure that the finite-size behavior of the density of contact points is 
indeed a consequence of conformal invariance, in Figure 11, we give the density of contact 
points Nj^\x) multiplied by [L sin(7ra;/L)/(l — cos(7rx/L))]i/'^ for different boundary rates 
a > 1, and a fixed, large value of L. One notices that within errors the value of f3 in eq. (l4.7l) 
stays unchanged. We repeated the study for a < 1 and for odd values of L. For L odd, as 
expected, we observe no a dependence of the density profiles. 
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Figure 10: (Color online) The density of contact points at x for a system of size L, Nj^ (x) 
multiplied by [Lsin(7rx/L)/(1 — cos(7rx/L))]-^/^ for a boundary rate a = 1 as a function of 
x/L for different lattice sizes L (even) = 512, 1024, 2048 and 4096. 



It is known that the finite-size scaling limit of the Hamiltonian spectrum which gives the 
time evolution of the system is independent of the boundary rate a [TB]. Consequently one 
could expect the space one-point function in the scaling limit to be also independent of a. 
This is indeed the case. 

To sum up, the density of contact points in the presence of a wall, has for any boundary 
rate a and for both L even and odd, the expression: 



N^l\x) - 0.8757 



1 - cos(^) 
Lsin(^) 



1/3 

(4.9) 



The same expression describes the density of contact points for the RPM {L odd) defined 
on one-step Dyck paths. 

From (14.91) one learns that for small values of x (large L) one has very few clusters. This 
can be explained by the existence of a large leftmost cluster (see Figure 8 and a quantitative 
argument below). At the other end of the system x = L — y, y small , the density of clusters 
is larger than in the Dyck path configurations (a factor of 0.8757 ■ 2^/'^ ~ 1.103 compared 
with 0.753). 

Up to now we have not discussed the exponent X = 1/3 and we are going to see something 
unexpected. According to (14. 4p one gets a scaling index Ap g = 1/6. This would imply that 
one has to chose p and q in equation (14. ip such that 

{3p-2qf = 5, (4.10) 

which is not possible for any rational values of p and q. On the other hand if one would 
take Ap^q = 1/3 and Api^gi = in eq. (l4.2p (both values compatible with (14. ip ) one would 
get X = 1/3 and s = 1/3, a non-zero value of the conformal spin and therefore a vanishing 
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Figure 11: (Color online) Density of contact points {x) multiplied by 

[(Lsin(7rx/L)/(1 - cos(7rx/L))] ^3 

function oi x/L for L (even) = 4096 and input 

rates a = 1, 2, 5, 10 and 20. 



density profile. We have not found an explanation for this puzzle. We are going to see that 
there is also a second one. 

2) Density of defects. 

We consider the RPMW in the defect configurations (see Appendix A and Figure 5). 
We have to stress that in the ballot path configurations, defects are nonlocal observables the 
same way as contact points are nonlocal observables in defect configurations. One should add 
that if one studies the XXZ quantum chain [B] , which has the same spectrum as the RPMW, 
both contact points and defects are nonlocal observables. In the case of the quantum chain 
the natural observable could be the local magnetization (erf) pT]. One should keep in mind 
that in the spin basis one looses the probabilistic interpretation of the ground state of the 
Hamiltonian. 

In the defect configurations, the RPMW Hamiltonian acts in the following way. The 
clusters evolve as in the RPM, adjacent defects {D) annihilate, defects hop over the adjacent 
clusters pealing them and moving them as a result of the hopping. What we have described 
up to now can be seen as a pair annihilation process 

D + D^^ (4.11) 

in an unquenched disordered media. This is the situation if the boundary rate a = and one 
has no defects in the stationary state. If the boundary rate is not zero, one injects defects 
in the system which compensate for the losses through annihilation and in the stationary 
states one finds defects and is interested in the density of defects profiles. 

For convenience, if a defect is at the point (x + 1/2), x an integer, we will shift it at x. 

We expect the density of defects D^^\x) to have an expression like (14.31) . As opposed to 
the density of contact points X is unknown, since the total number of defects, equal to the 
average height at the origin (ho) in the ballot path configurations, is not known. 
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Figure 12: (Color online) The density of defects D\\x) multiplied by Lsin(7rx/L) as a 
function of (l + cos(7rx/L))/2 for different lattice sizes. The profiles for even and odd sites 
are shown separately. 



We have observed that one has to consider separately the scaling properties of the density 
of defects on even and odd sites (denoted by d'^^'^{x) and D'f^^ix) respectively). Monte 
Carlo simulations on large lattices have shown that the scaling dimension of the local operator 
corresponding to the defects is: 

X = l, (4.12) 

and therefore one has: 

L sm(^) 

In Figurel2 one can see the scaling functions for an input rate a = 1. Surprisingly, G^^^° 
has a simple expression but not One observes that ^(^^^(1) = G(i)°(l) ~ 1.75. This 

implies that near the source we have the universal behavior: 

D^l\x) ~ ~ 0.577- (4.14) 

Tlx X 

for both even and odd sites, and that the total number of defects, equal to (/io)L^ in the 
ballot path picture, increases logarithmically with L. A careful counting gives: 

(/io)i^^ ~ 0.467 + 0.577 In L. (4.15) 

We have checked that the results presented above are independent of the boundary rate 
a as expected from conformal invariance. We now try to understand the exponent X = 1 
(see f lAj) ). This would imply (see (gT]) and flT2|) ): 

Ap,g = 1/2 and (3^-2^)^ = 13, (4.16) 

which as in the case of the density of contact points, is not possible. Similar to the previous 
case, one could take Ap ^ = 1, Ap/ g' = 0, both values compatible with (14.11) and get a nonzero 
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Figure 13: (Color online) The probability Fl{x) to have the first cluster ending at x in 
the RPM (L even) multiplied by x^/^, as a function of x/L. 



conformal spin. Notice that the 1/x fall off ( 14.14p of d'^1\x) is the mean-field result for the 
case in which the hopping is local and symmetric and the input of particles is also local [22] . 
The same behavior is obtained when fluctuations are taken into account [23]. In our case 
neither the source acts locally nor is the symmetric hopping local. 

3) Probability density for the first cluster. 

a) The RPM (L even). 

We consider the probability Fl{x) to have the first cluster end at the point x {i. e. the 
second contact point at x in a Dyck path). In Figure 13 we show how this function scales: 

FL{x)-f{f;)/x'/\ (4.17) 

This implies that the first cluster increases like 3.3 L^/^, similar to the average cluster size 
1/0.738488 - 1.354118^^/3 ^^^^ (g^j)). No surprises here. 

b) The RPM and the RPMW for L odd. 

This case is much more interesting. Firstly let us recall that F^\x) coincides with the 
probability of having the first cluster end at x in the one-step Dyck path configurations in 
the RPM model (L odd) (see fl3.10p and (13. lip ). The latter coincides with the of having 
the defect at (x + 1/2) in the one-defect configurations. This suggests that Fj^\x, L) might 
have other properties. 

Let us consider the following, normalized to 1, probability density function: 

Pi^^L) = Z..,i/3 ^ (4-18) 



L[sin(^' 



where 

C = 2V^r(|)/r(i) ~ 1.4936. (4.19) 
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Figure 14: (Color online) The probability distribution function to have the first cluster 
ending at x, for a system of size L, F^\x) multiplied by L(sin(7rx/-L)^/^ as a function of 
x/L for different values of L and boundary rates a. 



In Figure 14 we show that in the scahng Umit, Fi' (x) stays a probabihty density distribu- 
tion given by ( 14.18^ for any boundary rate a. This is unexpected and should be understood. 
Note the reappearance of the exponent 1/3 in (14.181) 

Notice that the first cluster is either small or very big {p{x, L) = p{L — x,L)) and that 
the average size of the first cluster is L/2 (this is a linear dependence as opposed to the L^/^ 
dependence in the RPM with L even). 

How one can derive a probability density (14.181) from conformal field theory is not clear 
to us. 

One of the most interesting aspects of our study is the deep connection between the 
RPM for L odd (one-step Dyck paths) and the RPMW model in the defect configurations 
(corresponding to ballot paths). The rightmost defect behaves like the single defect which 
defines the position of the leftmost cluster in the RPM (one-step Dyck paths picture). Like 
the single defect, the rightmost defect makes Levy flights and behaves like a "relativistic" 
random walker [5]. The physics to the right of the leftmost cluster is identical in both cases, 
therefore the density of contact points in the two models is the same. The difference between 
the two models is all in the first cluster. In the RPMW model one has defects in the first 
cluster, while in the RPM there are no defects. The content of the first cluster in the RPM 
can be easier understood in the one-defect picture in which the time-evolution operator acts 
in the left-right symmetric way. 

Before closing this section, let us discuss shortly the implication of the conjectures (13.161) 
for A^l\ the average number of local minima of the interface divided by L in the RPMW. 
In the large L limit one obtains 

'-^ 3/8 . (4.20) 

This number is identical with the one observed in the RPM (L even) [1]. This implies (see 
[1]) that the average number of tiles desorbed in avalanches is again 1.5. This coincidence 
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suggests that most of the desorption processes take place within the clusters. The different 
cluster structures in the RPM and the RPMW, do not affect most of the avalanches. This 
does not imply that the rare large avalanches [2] are the same in the two models. 

5 The raise and peel model with and without a wall in 
continuous time. Temperley-Lieb algebras 

In this section we show how the RPM and RPMW described in Section 2 were obtained. We 
also give the necessary background to derive the exact results given in Section 3 and derived 
in Sections 7 and 8. 

The continuous time evolution of a system composed of the states a = 1,2, ... ,n with 
probabilities Pa{t) is given by a master equation that can be interpreted as an imaginary 
time Schrodinger equation: 

j^p4t) = -j:H^Mt), (5.1) 

where the Hamiltonian H is an n x n intensity matrix: i?Q,/3 nonpositive {a ^ (3) and 
J2aHa,i3 = 0. —Ha,i3 is the rate for the transition \/3) \a). The ground-state eigenvector 
of the system |0), i^|0) = 0, gives the probabilities in the stationary state: 

\0)=J2Pa\a), P„=limP„(t). (5.2) 

a 

The normalization factor of the unnormalized probabilities Pa is (0| 0) where 

(0|=5:(a|, {0\H = 0. (5.3) 

a 

Since H is an intensity matrix, the real parts of the n eigenvalues Ea are nonnegative and, 
if complex, the eigenvalues come in conjugate pairs. The eigenvalue zero is not degenerate. 

For our present purposes we need the following observation ^: ii Qi {i = 1, . . . ,p) are 
generators of a semigroup algebra, then: 

H = j2a,{l-g,) (5.4) 

i=l 

acting from the left in the vector space defined by the regular representation or by a left 
ideal of the algebra, is an intensity matrix (oj are non-negative real numbers). 

In the following we are going to use the Temperley-Lieb (TL) and the one-boundary 
Temperley-Lieb semigroup algebras. 

The Temperley-Lieb algebra 7^ [9] is given by a set of (L—l) generators Ci, i = 1, . . . , L — 
1, subject to relations 

ef = Tei, ejCiiiCi = Ci, etej = ejCt if |2-j|>l. (5.5) 
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Here r is the parameter of the algebra which for our apphcations we put equal to one: r = 1. 
We notice that with this choice all structure constants of the algebra become units and the 
TL algebra becomes a semigroup algebra (see [3]). 

The one-boundary extension of the Temperley-Lieb algebra Tj^^^ (also known under the 
name of blob algebra [5]) is obtained by adding a boundary generator cq satisfying relations 

ejCo Vz > 1. (5.6) 



eiCoei = ei, 



We fix again the parameter uj = 1 and get a semigroup algebra. 

For our purposes it is convenient to represent the generators Cj and the boundary gener- 
ator Co as, respectively, tiles and a half-tile: 







eo 



(5.7) 

The tiles and the half-tile are placed on labeled vertical lines and they can move freely along 
the lines unless they meet other (half-)tiles: 



m 



1 2 3 ... i-1 i j+1 



1 2 3 ... i-1 i s+1 



Multiplication in the algebra corresponds to a simultaneous placement of several (half-)tiles 
on the same picture and an order of the product corresponds to moving the tiles downwards. 
Thus, the picture above corresponds to the commutation relations in (15.51) and (15. 6p . The 
remaining relations in (15.51) and (15. 6p are equivalent to the following pictures 





m 



(5. 



, . i— 1 i . . . 



. i—1 i ■ ■ ■ 



(5.9) 



1 2 ... 1 2 ... 



1 2 ... 12 ... 



We consider now the following projectors in the TL and the one-boundary TL algebras: 
Temperley-Lieb algebra and 



one-boundary TL algebra, 



L even: 



1 2 3 4 ... L-1 L 




one-boundary TL algebra L odd: 



606264 . . . Cl-I 



1 2 3 4 5 ... 1^1 L 
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We multiply the projectors by all elements of the respective algebras from the left (that 
is, dropping the (half-)tiles downwards) and generate in this way left ideals (left-invariant 
subspaces in the algebras). The projectors are non- invert ible and therefore the ideals are 
smaller then the algebras themselves. The dimensions of the resulting ideals are: Cl (see 
eq. (12. 2p ) for the TL algebra for L even, Cl+i for the TL algebra for L odd, and {^yL/2\) 
one-boundary TL algebra. These dimensions coincide with the numbers of the Dyck paths, 
one-step Dyck paths and the ballot paths, respectively. One can easily see a correspondence 
between the paths and the pictures of words in the ideals. The correspondence is illustrated 
below. 



Dyck path: 




646661636567 



One-step Dyck path: 





62616365626466 



1 2 3 4 5 6 L=7 



= 61606261636560626466 
12 3 ... L=T 1 2 3 4 5 6 L=7 

Obviously, the projectors correspond to the substrates in the Figures 1, 2 and 4 (see Sec. 2). 

The Hamiltonian governing the time fluctuations of the RPM and the RPMW are given 
by the following expressions (see (15.41) ) 

Hl = E(l - ^i) ' = «(1 - eo) + E(l - ei) , (5.10) 

i=l 1=1 

where a is the boundary rate. 

We now notice that the multiplication of the words from the ideals by the generators 6j 
and 60 drawn in pictures according to the rules (15.81) and (15.91) coincides with the elementary 
processes (adsorption, reflection and desorption) of the RPM and the RPMW models. 

The groundstate eigenvectors of the two intensity matrices (I5.10p have remarkable combi- 
natorial properties [HI [T21 [131 [IH [13 EH ES] • New combinatorial aspects of these groundstate 
eigenvectors are going to be shown in Sections 7 and 8 and will be used to derive the results 
enumerated in Section 3. 

One can use the representation of the generators 60, 6j in terms of the Pauli matrices 
to derive exact results. In this representation Hi and H^^^ are given by XXZ spin one-half 
quantum chain Hamiltonians with non-Hermitean boundary terms From the finite-size 
scaling limit of the Hamiltonian eigenspectra we can derive the conformal properties of the 
systems [TOl [6] . 



Ballot path: 
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6 Hexagon recurrence relations 



In this section we want to show that the expressions (13. ip and (13.21) for S{L, n) and P{L, n, m) 
are solutions of two recurrence relations. They can be seen as generalizations of the recur- 
rence relations obtained from the Pascal triangle. This section is pure mathematics. The 
connection with the RPM and RPMW is going to be made in Sections 7 and 8. 



6.1 Pascal's hexagon 

Let us place numbers on the sites of a trigonal lattice in such a way that any six numbers, 
occupying the vertices of an elementary hexagon as shown on the figure below, 



X 



a 



(6.i; 



satisfy the relation 

dc = a'a + h'h. (6.2) 

Introducing a coordinate frame {Ht , 1^} on the lattice (see figure above) and denoting S'(L, n) 
the number associated to a vertex (L"af + nlj) we can write the relation ( 16. 2p in the form 

S{L + l,n+l)S{L-l,n-l) = S{L - 1, n) S{L + 1, n) + S{L, n - 1) S{L, n + 1) . (6.3) 

This equation was obtained in p6] as a discretized version of the Boussinesq equation. It 
is a particular two-dimensional reduction of the three-dimensional discrete Hirota equation 
(see [271 EH] and references therein) which is also known as the octahedron recurrence in 
combinatorics (see [29] and references therein). 

We treat (16.31) as a recurrence relation which gives the number S{L+1, n + 1) (correspond- 
ing to c' in the figure) in terms of the numbers placed in the two upper rows (respectively, 
a, a', c and b') and at the left neighbour site from the same row (resp., b). We solve the 
recurrence relation by moving rightwards in a row and downwards row by row. Taking the 
following set of initial data 

5(L, -1) = 5(L,0) = 1, S{2n-l,n) = Vn > 0. (6.4) 

we can calculate numbers S{L, n) in a sector n> —1, L>2n — 1. As far as we know these 
boundary conditions were not considered before. 

The values of S{L,n) for small L and n are given in Figure 15. This figure and the 
recursive procedure described above reminds of the famous Pascal's triangle rule — an 
arrangement of the binomial coefficients in a triangle. A closer inspection of the solution 
adds more arguments in favour of such an analogy. First, we notice that all numbers on 
Figure 15 are integers, which is not trivial since solving the recurrence relation (16.31) one 
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1 ^ 1111111 1111 



n=0 

1 1111111 11111 




23456 789 10 11 

3 11 26 50 85 133 196 276 375 



26 170 646 1862 4508 9660 18900 

646 7429 45885 202860 720360 

45885 920460 9304650 

Figure 15: Pascal's hexagon rule. 

obtains ratios|§ Moreover, all the numbers in Figure 15 can be factorized into relatively 
small primes (compared to the numbers themselves). For example S{9,4), S'(10,5) and 
S'(12,5) are given respectively by: 

7429 = 17 ■ 19 ■ 23, 45885 = 3 ■ 5 ■ 7 ■ 19 ■ 23, 9304650 = 2 ■ 3^ ■ 5^ ■ 23 ■ 29 ■ 31. 

This factorization property suggests an idea to look for an expression of S{L,n) in terms of 
factorials, just as it is the case for the binomial coefficients. The solution of the recurrence 
relation (6.3), with initial condition (6.4), is given by (13. ip |14j . 

One can look for a more general solution of the Pascal's hexagon recurrence (16. 3p using 
the initial conditions depending now on a parameter a: 

^('^)(L,-1) = 1, S^''\L,0) = a^^-", S^''\2n-l,n) = Vn > 0. (6.5) 

Remarkably again, all the quantities S^^\L,n) in the sector n > —1, L > 2n — I, turn out 
to be polynomials in variable a with integer coefficients. These polynomials are particular 
examples of the polynomial solutions considered in the Appendix of Ref. [13] {S^°'\L,n) = 
FL,n{x = 1; 2/ = c^) ill the notation of [I3])- The polynomials S^""^ (L, n) can be also observed in 
the stationary states of the RPMW with an arbitrary boundary rate a (see Sec. 8). Obviously, 
one has S^^\L,n) = S{L,n). 

6.2 Split-hexagon 

We now introduce a vertical direction in the picture (16. ip and put numbers P{L, n, m) 
into correspondence with the vertices [Lit + nlj + rn^) of the 3(i lattice (see Figure 16). 

^ For a general discussion of such phenomena see Ref. [30] ■ 
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P(L-l.n,my. 



\ 5(L-l,n-l) i 
• .• • 



■ \S{L,n+l) 
• 



P(L,n-l,m-2) 



Figure 16: Split hexagon: the numbers S{L + l,n + 1), S{L — l,n) and S{L,n - 1) 
split, respectively, into the parts P{L + l,n + l,m), P{L — l,n, m) and P{L,n — l,m), 
m = 0, 1, . . . . 



We impose the following relations on P(L, ra, m)'s: 

n 

J2 P{L,n,m) = S{L,n), 

m=Q 

P{L,n,-k) = P{L,n,n + k) = > , n > , 

and the equations 

P{L + l,n + l,m) S{L-l,n-l) 

= P{L - 1, n, m) S{L + 1, n) + P(L, n - 1, m - 2) ^(L, n + 1) 



(6.6) 
(6.7) 



(6.^ 



Summing over m the relations (16.81) one obtains the hexagon relations (16.31) . therefore we 
call (16. 8p the split-hexagon relations. 

Taking into account (16.61) . the Eq. (16. Sp defines recurrence relations for the P{L,n,m)^s. 
Indeed, using (16.71) and the initial conditions (c.f. with the initial data (16.41) ) 



P(L,0,m) = 6, 



m,0 



P{L, -l,m) = 5, 



m,—l ) 



and P{2n,n,m) = Vm > 1, (6.9) 



one can calculate the numbers P{L,n,m) starting from the row n = 1 downwards and 
moving rightwards from the point L = 2n in a row. For each given pair L and n one 
calculates P{L, n, m) for all m = 0, 1, . . . n. 

We notice that the split-hexagon recursion is a kind of "bootstrap" procedure as after 
each step (for example, the calculation of P{L, n + 1, m) for all m = 0, 1, . . . , n + 1) one has 
to find a coefficient for the equations which are solved at the next step of the recursion (that 
is, the coefficient S{L,n + 1) = Y1^=q P{L, n + l,m) in the equations for P{L + 1, n + 1, m)). 
We also remark that there is no need to set the boundary values 

P{2n,n,0) = S{2n,n-1) and P{2n,n,l) = 

as they are completely defined by the data from the previous (n — l)-th row only (the last 
term in the split-hexagon relation (16. 8p vanishes for m = or 1). 
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The relations (16.81) are new and can be generalized changing the condition (16. 6p . The 
continuous version of fl6.8p is not yet known. 

The solution of the split-hexagon recursion is given by integer numbers with simple 
factorization properties, e.g., 

P(12,4,2) = 2041020 = 22-32-5- 17-23-29 , P(12,4,3) = 2281140 = 22-32-5- 19-23-29 . 

This suggests the P(L,n, m)'s can be written as a product of factorials. One can show that 
the expression given by (13. 2p for P[L,n,m) is the solution of the split-hexagon recurrence 
relations (16. 6p and (16. 8p with the boundary conditions (16. 7p and (16. 9p . 

7 Combinatorics in the stationary states of the RPM 

We have shown in the last section that the quantities P{L, n, m) are solutions of the split- 
hexagon recurrence relations (16. 8p and we gave their expression (13. 2p . We now make the 
connection with the stationary states of the RPmH 

The stationary state of the RPM for a system of size L is described by the eigenvector 
|0)j^ of the operator Hl (15.100 corresponding to eigenvalue zero. It is a linear combination 
of the words in the Dyck path representation (see Sec. 5) defined up to a common constant 
factor. A suitable normalization for |0)^ is to set the coefficient of the maximal Dyck path 
(see figures below) to be equal to 1. 



Maximal Dyck paths: 




for L odd. 



1 L=6 1 L=5 



This is the smallest coefficient in the vector |0)^ and the corresponding profile is the least 
probable in the stationary state. With this choice, all the weights of the different Dyck paths 
in the ground-state eigenvectors are integer numbers which can be related to an enumeration 
of loops in a fully packed loops ice- model [24]. The most probable stationary configuration 
is the substrate shown as a dashed region on the figures abovefl 

Consider a Dyck path with L sites. Let /i be a number from the list 

0, 1,..., L^J, 

such that the Dyck path has no local minima (return points) below the height h. (See 
Figure 17). Let us denote by m (m > 0) the number of returns at the height h for a given 
path. We consider all the configurations which have m return points at the height h and the 
remaining ones higher. 

^ Actually it was by observing the combinatorial properties of the RPM's stationary states that we were 
led to the hexagon and the split-hexagon relations. 

^Intuitively, the more reflection points (local maxima) one has in a profile, the larger is its contribution 
to the stationary state. 
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L=16 



Figure 17: (Color online) Examples of Dyck paths of the size L = 16 satisfying the 
selection criteria for /i = 2 n = 5). The number m of return points at the height 
h = 2 (red points on the pictures) is equal, respectively, to 5, 2, and 0. 

Conjecture. The sum of all the coefficients in |0)^ corresponding to these configurations 
is given by P{L,n,m) . Here 

n = [^^^y^J - /i , = 0, 1, . . . , [^^y^J , and m = 0, 1, . . . , . (7.1) 

n counts the height from the top of the triangle instead of its basis. 

This conjecture was checked for lattices of size L < 13. There are more consistency 
checks for this conjecture: 

• One clearly has P{L,n + 1,0) = J2m=o ^i^^^y''^) — S{L,n) (see (16.61) ). which is 
the sum of all the coefficients in |0)^ corresponding to the configurations with local 
minima at values higher or equal to h = {\_{L — 1)/2J — n) and is given by (13. ip . 
This expression was conjectured earlier and checked independently on small lattices in 
Refs. [131 HU. 

• More importantly, it was proven analytically in Refs. [311 [32] that S{L) (see (13. 3p and 
(13.41) ) coincides with the normalization factor (sum of all the weights) of |0)^: 

^(0|0)^ = S{L) = S{L,[{L -l)/2\). (7.2) 

The expression for the probability of having k clusters for a system of size L given in 
Sec. 3 (see (13.51) ) is an obvious corollary of the conjecture. 

8 Clusters in the stationary states of the RPMW and 
RPM 

In this section we prove the expressions (13. 6p and (13.70 for the probabilities to have k clusters 
for a system of size L as well the relations (I3.10p - (l3.13p for the probability densities observed 
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in the RPMW. To this end we make use of the algebraic background (the TL algebras 
described in Section 6), and of the conjecture presented in Section 7. 

The stationary state of the RPMW for a system of size L is described by the ground 
state eigenvector of the operator H^^^ (15.101) which is a linear combination of the words in the 
ballot path representation (see Sec. 5) with coefficients (weights) depending on the boundary 
rate a. We use a notation lO)^-* for this eigenvector and, if necessary, we explicitly show the 
value of a in the notation, say |0)^^^ for a = 1. A suitable normalization for |0)^'*^ is to set 
the coefficient of the maximal ballot path (see figure below) to be equal to a^^^^J. In case 
a = 1 it is the smallest coefficient in |0)^^'*. 



Maximal ballot path: 




With this choice, all the weights of the different ballot paths in the ground-state eigenvectors 
become polynomials in a with integer coefficients which can be related to a "weighted" 
enumeration of loops in a fully packed loops ice- model [I5] . As it was proven in [32] , in the 
case a = 1, the normalization factor (sum of all weights) of |0)^^^ is equal to the number of 
vertically and horizontally symmetric alternating sign matrices of a size (2L + 3) (see (13.31) 
and m 

^'M 0)i'^ = ^2^1+3 = S{L) S{L + 1) . (8.2) 
For an arbitrary boundary rate a, the normalization factor reads (see [32], Eq. (4.8)) 

('^2(01 0)i"^ = S{L) S^''\L + 1) , (8.3) 

Here S^^^L) = S^^^L, [(L - 1)/2J), and S^''\L,n) is a polynomial solution of the Pascal's 
hexagon recurrence (16.31) with the initial conditions (16.51) . 

Now we are going to use homomorphisms from the one-boundary TL algebra T^^^^ to the 
TL algebras 7^ and T^+i to establish relations between the ground state vectors of the RPM 
models with and without a wall. One has to consider the cases L even and odd separately. 

8.1 Case of L odd. 

We consider the homomorphism between the one-boundary TL algebra and the TL algebra 
(see (15.51) and (15. 6p ) of the same size L, defined by 

Ti'""^ ^ Tl: eo 1, ^ e, = 1, . . . , L - 1. (8.4) 

This homomorphism maps the intensity matrix of the RPMW (a is arbitrary) onto the 
intensity matrix of the RPM for the same number of sites L (see (I5.1UI) ) 

^ Hl , (8.5) 
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Figure 18: Four different ballot paths (shown on top) map to the same one-step 
Dyck path (shown below) under the homomorphism (18.41) . Marked by asterisk "*" 
are the half-tiles which one erases during the mapping. Marked by dagger "f" are 
the tiles which one then removes using the rule (15.81) . There are r = 2 reflection 
points at height h = 1 in the first cluster of the one-step Dyck path (see marked 
points). Therefore, one finds 2^ = A different ballot paths in the preimage. The 
starting points for the layers of tiles are marked on the pictures of ballot paths. 



and so should map their unique groundstate eigenvectors: 

|0)ir) ^ ^('^)(L + 1)|0),, (8.6) 

where we have computed the scaling constant S^^-^L + 1) by comparing the expressions for 
normalization factors (17.21) and (18. 3p . 

Let us consider in more details an action of the map (18.61) on a ballot path. Graphically 
one can realize this mapping using the following procedure (see Figure 18): one first erases 
all the half-tiles from the picture of the ballot path and then, one reduces the resulting 
expression using the rules (15.81) (actually one needs only the last equality from (15. 8p ). There 
can be several ballot paths mapping onto the same one-step Dyck path under this procedure. 
One can reconstruct the corresponding ballot paths by dressing the one-step Dyck path with 
the layers of tiles starting at any of its return point (local minima) at height h = 1 in the first 
cluster and ending by a half-tile at the wall (see Figure 18). Assuming there are r such return 
points in the one-step Dyck path, one can construct 2*" different "dressed" ballot paths. 

We observe that the profiles of all the ballot paths mapping under (18. 4p onto the same 
Dyck path coincide except for their leftmost clusters. They are identical to the profile of the 
corresponding Dyck path, again, with an exception of its leftmost cluster. We now notice 
that, according to (18. 6p . in the RPMW (L odd) the stationary weights of the ballot paths 
mapping to the same Dyck path sum up (modulo the constant factor) S^'^^L + 1) to the 
stationary weight of the corresponding Dyck path in the RPM (same size L). These two 
observations justify the following statement. 

Proposition 1. Consider the RPM and the RPMW for the same number of sites L (L odd) 
and any boundary rate a. For both models the stationary probabilities to have the first cluster 
at distance x and any given configuration to the right of it are the same. 
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L=10 



i=ll 



Figure 19: A ballot path which has height at the wall maps under the homomor- 
phism (18.71) to the one-step Dyck path with hi = 0. There are no other ballot paths 
mapping to this Dyck path. 



In particular, for L odd, the relations (I3.10p and (13.121) are valid and the conjecture 
announced on page [26] implies the equality (13.61) . 

8.2 Case of L even. 

In this case we restrict our consideration to the case a = 1 only. We consider another 
homomorphism mapping the one-boundary TL algebra of size L onto the TL algebra of size 
(L + 1): 

j-(iB) _^ j-^^^ . ^ ^ ^^^^ Vi = 1, . . . , L - 1. (8.7) 

This homomorphism maps the intensity matrix of the RPMW, a = 1, to the intensity matrix 
of the RPM 

H^^ ^ Hl+1 , (8.8) 
and so should map their unique ground-state eigenvectors: 

|0)W ^ S{L) |0),^, , (8.9) 

where we have computed the scaling constant S{L) by comparing the expressions for their 
normalization factors (17.21) and (18. 2p . 

With respect to the action of homomorphism (18.71) the set of even size L ballot paths 
separates naturally into two parts. 

a) The ballot paths which have zero height at the origin, = 0, map onto the one-step 
Dyck paths with hi = (see Figure 19). The mapping (18.71) on these subsets of paths is a 
one to one correspondence. 

b) The ballot paths which have heights /iq = 2, 4, . . . , ^ at the origin map onto the one-step 
Dyck paths with hi = 2 (see Figure 20). Graphically one can realize the map (18. 7p for such 
paths in the following way: firstly, one completes each half-tile on the picture of the ballot 
path to a full tile and then, one reduces the resulting expression using the rules (15.81) (one 
needs only the last equality from (15.81) ). This map is not a one to one correspondence. The 
number of ballot paths in the preimage of a given restricted, hi = 2, one-step Dyck path is 
equal to 2^, where r is the number of return points at height /i = 1 in the first cluster of 
the one-step Dyck path, the point h^ = 1 is excluded (see Figure 20). The reconstruction of 
the preimage in this case is obtained again, by dressing the one-step Dyck path with layers 
touching the wall (see caption of Figure 20). 

We observe that in both ) and b), the profiles of the ballot paths, except their 

leftmost clusters, are identical to the profiles of the corresponding one-step Dyck paths. We 
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12 3 L=ll 



Figure 20: Four different ballot paths (shown on top) map to the same one-step 
Dyck path (shown below) under the homomorphism 08. 71) . The half-tiles on the 
pictures of ballot paths are completed to full tiles. Marked by dagger " f" are the 
tiles which one removes during the mapping using the rule (15.81) . There are r = 2 
reflection points at height h = 1 {ho excluded) in the first cluster of the one-step 
Dyck path (see marked points). Therefore, one finds 2^" = 4 different ballot paths in 
the preimage. The starting points for the layers of tiles are marked on the pictures 
of ballot paths. 



also notice that the leftmost clusters become one step larger under the map (18. 7p . Therefore, 
taking into account the relation (18. 9p we conclude: 

Proposition 2. For the RPMW of even size L and the boundary rate a = 1, the stationary 
probability of having the first cluster at distance x and any given configuration to the right 
of it is equal to the stationary probability for the RPM of odd size (L + 1), to have the first 
cluster at distance [x + 1) and the same configuration to the right of it. 

In particular, the relations (13. lip and (13.131) are satisfied. 

Finally, we demonstrate that the proposition 2, together with the conjecture presented in 
page [26] implies the equality (13.70 . Taking into account that the ballot paths with k clusters 
project under (18. 7p to the one-step Dyck paths with ik + l) clusters in case a) (see Figure 19) 
and with k clusters in case b) (see Figure 20) one can write 

L even: P^j^\k) = Pl+i(A;, hi = 2) + PL+i{k + 1, = 0) 

= Pl+i(A:) - Pl+i(A;,/ii = 0) + Pi+i(A; + l,/ii =0), (8.10) 

where PL+i{k, hi = 0(2)) are stationary probabilities to find k clusters and the height hi = 0(2) 
in the first site for the RPM with the odd number of sites (L + 1): 

PL+iik, hi = 0) + PL+i{k, hi = 2) = PL+i{k) . 

The first term in the sum (I8.10p is a stationary contribution of the profiles with {k — 1) 
return points at height h = (<(=> n = L/2, see Figure 17) which is given by the conjecture: 

PL+i(k) = 1 P(L + l,f,fc-l). (8.11) 
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To calculate the last two terms in the sum flS.lOl) we take into account the left-right symmetry 
of the one-defect profiles in the "defects" picture of the one-step Dyck paths (see Appendix 
A). As it is seen on the following picture 



One-step 
Dyck paths: 




Single defect 
configurations: 





in the RPM of the odd size (L -|- 1) the stationary weights of profiles with {k — 1) contact 
points [i.e. with k clusters), one of them is at site 1 (/ii = 0), are equal to the stationary 
weights of profiles with no contact points, but with (fc — 2) return points at height h = 1 
(■v^ n = "I — 1, see Figure 17). Therefore 

PL+i{k, /ii = 0) = P(L + 1, f - 1, A; - 2) . (8.12) 

From fl8.10p - fl8.12p we obtain for L even 

PL\k) = ^(^{P(L + l,f,A;-l)-P(L + l,f -l,A;-2) + P(L + l,f -l,fc-l)}, 

which is equivalent to (13.70 provided one takes into account the identity 

P(L + l,f,m)-P(L + l,-|-l,m-l) = P(L + 1, f , m + 1) for L even and Vm = 0, . . . , f . 



9 Conclusions 

As far as we know the RPM with and without a wall is the only example of a one-dimensional 
fluctuating interface which is conformal invariant (the central charge of the Virasoro algebra 
c = since in a stochastic process the ground-state energy is zero for any system size). 
The time dependence of average quantities are under control since the finite-size scaling 
Hamiltonian spectrum is known (see Refs. [ini[16]). For space dependent phenomena, the 
situation is more complex since one doesn't know which critical exponents appear. What 
we do know is that, in the finite-size scaling limit, the functional dependence of density 
profiles in the stationary state are fixed once the scaling dimensions of the local operators 
and the boundary conditions are specified (see (14. 3 p and (14.50 ). Our aim was to identify 
local operators by checking if the expressions (14.31) and (14. 5 p are confirmed by the data. 

First we have considered the contact-point density profiles in the RPM and in the RPMW 
and found the density profiles in the finite-size scaling limit using Monte Carlo simulations. 
They have the expected behavior (see (14.30 . (14.60 and (14.100 ) for an operator with scaling 
dimension X = 1/3. This value is puzzling since as explained in Section 4, it is incompatible 
with the Kac table for an operator of conformal spin zero (see fl4.10p ). The material presented 
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in Sections 6, 7 and 8 in which one computes the probabihties to have k clusters {k+1 contact 
points) in the system gives a rigorous proof that the exponent is indeed X = 1/3. We would 
like to qualify what we mean by rigorous. The derivation is based on finding solutions of 
the split-hexagon recurrence relations (16.81) for the quantities P{L^m^n). The identification 
of these quantities with certain Dyck path configurations as explained in Sec. 7 remains a 
conjecture. 

We have also considered the defect density profiles in the RPMW model and found the 
expected finite-size scaling function with an exponent X = 1 which is also puzzling since 
it is also incompatible with the Kac table for a spinless operator (see (14.161) ). We have 
no rigorous proof that the exponent is 1 but the physics which follows (see (14.141) ) is very 
plausible. 

One easy way out of these puzzles is to bring the argument that in a c = theory, 
everything is possible. We are tempted by another scenario: in both cases one has the 
conformal spin s = X. In this case both values 1/3 and 1 are given by the Kac table but we 
run into a different problem: an operator with conformal spin has a vanishing density. Why 
this does not happen in our case, is another puzzle. Possibly there is a simple explanation, 
but it escaped us. 

The model has many interesting properties as discussed in detail in Section 4. One of 
them looks to us really unexpected: the connection between the clusters in the RPM and 
the parameter dependent (the boundary rate a) RPMW for an odd number of sites. Except 
for the first cluster, the two models are identical. The proof of this statement is in Section 
8. It is perhaps no accident that the function which gives the probability of finding the first 
cluster at a distance x from the wall has special properties (see (I4.18P ). The implications of 
these observations on time dependent average quantities will be discussed elsewhere. 

We would like to comment about the content of the mathematical part of this paper. Its 
relevance goes beyond finding a critical exponent and mapping of two models. In Section 6 we 
give solutions for the so called Pascal's hexagon bilinear recurrence relations corresponding 
to certain boundary conditions. We also introduce new bilinear split-hexagon recurrence 
relations and solve them for certain boundary conditions. Their connection to integrable 
nonlinear differential equations is still to be found. 
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A Appendix. 

The RPM as a pair annihilation process with a source 

In Sec. 2 we have shown that the one-step configurations can be mapped onto configurations 
with one defect and clusters (see Figures 2 and 3) and that the ballot path configurations 
can be mapped onto defect configurations. In a defect configuration, clusters start and end 
on sites and defects occupy links. We give some simple examples of the mapping. 

In Figure 21 we show the three ballot paths for L = 3 and in Figure 22 we show the 
corresponding three defect configurations. The mapping is shown in Figures 23 and 24 for 
L = 4. For a defect configuration with L even (odd), one has an even (odd) number of 




12 3 



Figure 21: The three ballot path configurations for L = 3 (four sites). In the stationary 
state of the RPMW all configurations are present while in the case of the RPM only 
configurations (a) and (b) (one-step Dyck paths) occur. The configurations should be 
compared with these of Figure 22 



defects. 

The dynamics of the RPMW, given in Sec. 2 describing the evolution of the interface 
defined by ballot paths, can be translated in the evolution of the system defined by defect 
configurations. The evolution of the clusters is the same as in the RPM but the defects 
(impurities) hop over the adjacent clusters and in the hopping process they peel the clusters. 
When two adjacent defects annihilate they are replaced by a small cluster "building" the 
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1 2 3 1 2 3 




Figure 22: (Color online) The three configurations for L = 3 (four sites) formed by 
defects (red circles) and RSOS clusters. The corresponding ballot path configurations in 
the RPMW are given in Figure 21. 




01234 01234 




01234 01234 




1 2 3 4 1 2 3 4 



Figure 23: The six ballot path configurations for L = 4 (five sites). In the stationary 
state of the RPMW all configurations are present while in the case of the RPM only 
configurations (a) and (b) occur. The configurations should be compared with these of 
Figure 24 



substrate. The source acts nonlocally: it adds a defect on the first site and another one at 
the end of the first cluster. We show now in detail the evolution process. 

We consider a rarefied gas of tiles falling on a defect configuration with L + 1 sites. 
With a probability pb = 1/{L — 1 + a) a tile hits the bulk of the configuration (the sites 
z = 1, . . . , L — 1). With a probability ps = a/ {L — 1 + a) a. tile hits the boundary site i = 0. 
a is called the boundary rate. 

The evolution rules depend on where the tile hits the configuration. In Figure 25 we 
show the different cases. In the cases a, b and c the tile hits a cluster and the rules are the 
same as in the RPM given in Sec. 2. We give the rules for the cases d and e (see Figure 25) 

(d) The tile hits the site i which is at the right end of a cluster hj > hi^c = hi = 
(j = i — c + 1, . . . , 2 — 1). A defect is present on the link {i, i + 1). The defect hops 
on the first left link of the cluster {i — c,i — c + 1) . In order to create a free link 
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(c) 



(e) 




(d) 



(f) 



12 3 4 



12 3 4 



Figure 24: (Color online) The six configurations for L = 4 (five sites) formed by defects 
(red circles) and RSOS clusters. The corresponding ballot path configurations in the 
RPMW are given in Figure 23. 



Figure 25: (Color online) One of the defect configurations for L = 21 (22 sites). There are 
3 defects (red circles on the links) and 3 clusters. Also shown are six tiles (tilted squares) 
a - f belonging to the gas. When a tile hits the surface the effect is different in the six 
cases. 



on which the defect hops, part of the outer layer is desorbed: if hj^i — hj is negative 
{i — c+l<j<i — 2) the tile, with the center at coordinates (j, hj — 1), evaporates 
(see Figure 26a), afterwards, the remaining cluster moves to the right by one lattice 
spacing {hj — /ij+i), as in Figure 26b. An alternative way to visualize the process 
is shown in Figure 26b. First a whole layer of the cluster evaporates hj hj — 1 
[j = i — c+1, . . . — The defect hops from the link [i, i + to the link {i — c, i — c+1) 
and a half-tile cluster (part of the substrate) is created: /ij-i = 0,hi = 1, /ij+i = 0. 

(e) The tile hits the site i between two adjacent defects placed on the links {i — 
and + 1). The two defects annihilate and in their place appears a small cluster 
{hi_i = hi^i = 0,hi = 1) (see Figure 27). 

Up to now we have described the action of the RPM Hamiltonian (see (5.10)) on the 
defect configurations. This model was used in [ARl] to study the pair annihilation process 




12 3 



21 



35 



12 3 



21 




12 3 



Figure 26: (Color online) The new profile after the tile d in Figure 25 has hit the surface 
at the right end of a cluster. The defect hops to the left end of the cluster, first peeling a 
layer of 5 tiles (a), and next, the peeled cluster is translated by a lattice spacing unit to 
the right (b). 




12 3 ^ - 21 

Figure 27: (Color online) The new profile after the tile e has hit the surface between two 
defects. The defects have disappeared and in their place one gets a new small cluster. 

of defects in an unquenched disordered system. We now consider the effect of the boundary 
operator cq in the Hamiltonian H^^^ (see (5.10)). This corresponds to the effect of a tile 
hitting the boundary site i = (case / in Figure 25). 

(f) If the tile hits the site i = and the link (0, 1) is occupied by a defect, the tile is 
reflected. If there is no defect on the first link (this is the case / in Figure 25), the 
site i = is the left-end of a cluster hj > = he = {j = 1 , . . . , c — 1) two defects 
are added to the configuration (one defect on the link (0, 1) another one on the link 
(c — l,c)) after peeling and moving the first cluster (see Figure 28a). If hj+i — hj is 
positive (0 < j < c — 3) the outer tile, with center at coordinates (j + 1, hj^i — 1), 
evaporates, and the remaining cluster is moved to the left by one unit. An alternative 
way to visualize the effect of the tile which acts as a source of defects, is to first desorb 
a layer from the first cluster hj — > hj — 1 {j = 1, . . . , c — 1). Next the defects are 
introduced on the two links which became free as the result of the desorption process 
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(see Figure 28b). 




12 3 21 



Figure 28: (Color online) The new profile obtained after a tile / in Figure 25 has hit the 
site i = (the left side of a cluster defined by two contact points Hq = he = 0, c = 12 
in Figure 25). One first peels a layer of the cluster (5 tiles in the figure) (a), moves the 
remaining cluster by one unit lattice spacing to the left and then adds two defects (blue 
circles) (b). One on the link (0, 1) and one on the link (c — 1, c). 



References 

[1] J. de Gier, B. Nienhuis, P.A. Pearce and V. Rittenberg, 'The raise and peel model of 
a fluctuating interface', J. Stat. Phys. 114 (2004) 1-35 ^ arXiv: cond-inat/0 301430i 

[2] F.C Alcaraz, E. Levine and V. Rittenberg, 'Conformal invariance and its break- 
ing in a stochastic model of a fluctuating interface'. J. Stat. Mech. (2006) P08003 



arXiv : cond-mat/0604223 , 



[3] F.C. Alcaraz, and V. Rittenberg, 'Different facets of the raise and peel model'. 
J. Stat. Mech. (2007) P07009 [arXiv : cond-mat70703725| 

[4] J. de Gier, 'Loops, matchings and alternating- sign matrices', Discr. Math. 298 (2005) 
365-388 larXiv : math . CO/021 1285] 

[5] P.P. Martin and H. Saleur, Lett. Math. Phys. 30 (1994) 189-206, 'The blob algebra and 
the periodic Temperley-Lieb algebra', arXiv:hep-th/9302094 , 

[6] A. Nichols, V. Rittenberg, and J. de Gier, 'One-boundary Temperley-Lieb algebras in 
the XXZ and loop models', J. Stat. Mech. (2005) P05003 arXiv: cond-mat/0411512, 



37 



[7] K. Shelton, 'The Singled Out Game', Mathematics Magazine, 78 (2005) 15. 

[8] F.C.Alcaraz and V.Rittenberg, 'The pair annihilation reaction D + D ^ in 
disordered media and conformal invariance', Phys. Rev. E 75 (2007) 051110 
arXiv : cond-mat/0612272l 

[9] H.N.V. Temperley and E.H. Lieb, 'Relations between percolation and colouring prob- 
lems and other graph theoretical problems associated with regular planar lattices: some 
exact results for the percolation problem', Proc. Roy. Soc. London Ser. A 322 (1971) 251- 
280. 

[10] H. Saleur and M. Bauer, 'On some relations between local height probabilities and 
conformal invariance', Nucl. Phys. B320 (1989) 591-624. 

[11] A.V. Razumov and Yu.G. Stroganov, 'Spin chains and combinatories', 
J. Phys. A: Math. Gen. 34 (2001) 3185-3190 | arXiv:math. CO/0012141 1 

[12] A.V. Razumov and Yu.G. Stroganov, 'Combinatorial nature of ground state vector of 
0(1) loop model, Theor. Math. Phys. 138 (2004) 333-337 arXiv:math. CO/0104216. 

[13] S. Mitra, B. Nienhuis, J. de Gier and M.T. Batchelor, 'Exact expressions for correla- 
tions in the ground state of the dense 0(1) loop model'. J. Stat. Mech. (2004) P09010 
arXiv : cond -mat/0401245 , 

[14] P. Pyatov, 'Raise and Peel Models of fluctuating interfaces and combinatorics of Pas- 
cal's hexagon'. J. Stat. Mech. (2004) P09003 arXiv : math-ph/0406025, 

[15] J. de Gier, V. Rittenberg, 'Refined Razumov- Stroganov conjectures for open bound- 
aries', J. Stat. Mech. (2004) P09009 .arXiv : fflath-ph/0408042[ 

[16] J. de Gier, A. Nichols, P. Pyatov and V. Rittenberg, 'Magic in the spectra of the XXZ 
quantum chain with boundaries at A = and A = —1/2', Nucl. Phys. B729 (2005) 
387-418 a rXiv : hep-th/0505062, 

[17] D.M. Bressoud, 'Proofs and Confirmations. The Story of the Alternating Sign Matrix 
Conjecture', 1999 Cambridge University Press, Cambridge 

[18] G. KvLpeiheig,' Symmetry classes of alternating- sign-matrices under one roof, 
Ann. of Math. 156, no. 3 (2002) 835-866 [arXiv : math . C O/00081841 

[19] D.P. Robbins, 'Symmetry classes of alternating sign matrices, 2000 
larXiv : math . CO/0008045 [ 

[20] T.W. Burkhardt and T. Xue, 'Density profiles in confined critical systems and con- 
formal invariance', Phys. Rev. Lett. 66 (1991) 895-898; Nucl. Phys. B354 (1991) 
653-665. 

[21] I. Affleck, 'Edge magnetic field in the xxz spin-\ chain', J. Phys. A: Math. Gen. 31 
(1998) 2761-2766. 



38 



Z. Cheng, S. Redner and F. Leyvraz, 'Coagulation with a steady point monomer source', 
Phys. Rev. Lett. 62 (1989) 2321-2324. 

H. Hinrichsen, V. Rittenberg and H. Simon, 'Universality properties of the stationary 
states in the one- dimensional coagulation- diffusion model with external particle input', 
J. Stat. Phys. 86 (1997) 1203 larXiv : cond-iiiat/9606088| 

M.T. Batchelor, J. de Gier, B. Nienhuis, 'The quantum symmetric XXZ chain at A = 
— 1/2, alternating sign matrices and plane partitions', J. Phys. A: Math. Gen. 34 
(2001) L265-L270 [arXiv : cond-mat/0101385| 

P.A. Pearce, V. Rittenberg, J. de Gier, B. Nienhuis, 'Temperley-Lieb Stochastic Pro- 



cesses', J. Phys. A: Math. Gen. 35 (2002) L661-L668 arXiv:math-ph/0209017 



E. Date, M. Jimbo, and T. Miwa, 'Method for generating discrete soliton equations III, 
Journ. Phys. Soc. Japan 52, No.2 (1983) 388-393. 

R. Hirota, 'Discrete analogue of a generalized Toda equation', Journ. Phys. Soc. Japan 
50 (1981) 3785-3791. 

A. Zabrodin, 'A survey of Hirota 's difference equations', Theor. Mat. Fiz. 113 (1997) 



1347-1392 larXiv : solv-int/9704001t 



D.E. Speyer, 'Perfect Hatchings and the Octahedron Recurrence', Journal of Alge- 
braic Combinatorics, 25, no.3 (2007) 309-348 [arXiv : math . CO/04024521 



S. Fomin and A. Zelevinsky, 'The Laurent phenomenon'. Advances in Applied Mathe- 
matics 28, no.2 (2002) 1 19-144 jarXiv: math rC0/010424Tl 

P. Di Francesco, 'Inhomogeneous loop models with open boundaries'. 



J. Phys. A: Math. Gen. 38 (2005) 6091-6120 arXiv : math-ph/0504032 



P. Zinn- Justin, 'Loop model with mixed boundary conditions, qKZ equation and alter- 
nating sign matrices', J. Stat. Mech. (2007) P01007 iarXiv:math-ph/0610067. 



39 



